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ABSTRACT 

The fraction of compact active galactic nuclei (AGNs) that exhibit interstellar scin- 
tillation (ISS) at radio wavelengths, as well as their scintillation amplitudes, have been 
found to decrease significantly for sources at redshifts z > 2. This can be attributed to 
an increase in the angular sizes of the /ias-scale cores or a decrease in the flux densities 
of the compact ^as cores relative to that of the mas-scale components with increasing 
redshift, possibly arising from (1) the space-time curvature of an expanding Universe, 
(2) AGN evolution, (3) source selection biases, (4) scatter broadening in the ionized 

■ intergalactic medium (IGM) and intervening galaxies, or (5) gravitational lensing. We 

■ examine the frequency scaling of this redshift dependence of ISS to determine its origin, 
using data from a dual-frequency survey of ISS of 128 sources at < z < 4. We present 

Qh! a novel method of analysis which accounts for selection effects in the source sample. 

We determine that the redshift dependence of ISS is partially linked to the steepen- 
+j3 ; ing of source spectral indices (af g) with redshift, caused either by selection biases or 

AGN evolution, coupled with weaker ISS in the af; 4 , < —0.4 sources. Selecting only the 
—0.4 < af; 4 , < 0.4 sources, we find that the redshift dependence of ISS is still significant, 
but is not significantly steeper than the expected (1 + z) 5 scaling of source angular sizes 
£T) ■ due to cosmological expansion for a brightness temperature and flux-limited sample of 

sources. We find no significant evidence for scatter broadening in the IGM, ruling it out 
| as the main cause of the redshift dependence of ISS. We obtain an upper limit to IGM 

\q | scatter broadening of < 110//as at 4.9 GHz with 99% confidence for all lines of sight, 

and as low as < 8uas for sight-lines to the most compact, ~ lOuas sources. 
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Introduction 



Radio waves originating from galactic and extragalactic sources are scattered as they propagate 
through the ionized interstellar medium (ISM) of our own Galaxy, distorting the wavefronts and 
generating interference patterns on the plane of the Earth. Relative motion between the Earth 
and the scattering region causes these interference patterns to drift across the observing telescope, 
and are detected as flux density variations in the source. This phenomenon, known as interstellar 
scintillation (ISS), has been observed to modulate the flux densities of very compact radio sources 
such as pulsars and quasars, and has been used as a pow erful probe of th e phys ics of the ISM and 
of the background sources themselves (see iRickettl (|1990l ) and lLazio et al.l (|20Q4j ) for reviews) . 



The presence of ISS is highly sensitive to the angular sizes of the sources. In an extended 
source, its amplitude is suppressed relative to that of a point source. This is because the scattered 
wavefronts at the observer's plane smear out when integrated over the emission originating from 
each element of the extended source, similar to the quenching of atmospheric scintillation at optical 
wavelengths in planets whose angular sizes are larger than that of stars. Since the interfering waves 
arrive from regions separated by distances up to an order of ~ 10 6 km apart at the scattering 
cloud, the ISM essentially functions as an interstellar interferometer, allowing ISS to probe source 
sizes down to //as scales. This provides a resolution orders of magnitude better than ground-based 
VLBI. In fact, its counterpart, interplanetary scintillation (IPS), has been used in the past to 
deter mine the angular sizes of radio sources at sub-arcsecond scales prior to the development of 
VLBI kittle fc Hewishlll96fil : Icohen et alJll967h . 



In the last few decades, the body of evidence linking the intraday variability (IDV) observed 
in many compact, flat-spectrum Active Galactic Nuclei (AGNs) at cm wavelengths to the phys- 
ical process of ISS has grown considerably. T ime delays between v ariability patterns have been 
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2003), caused by changes in 



orbital motion of the Ear th. The 4.9 GHz, four-epoch Micro-arcsecond Scintillation-Induced Vari- 



ability (MASIV) Survey (jLovell et al.l 120031 ) found 58% of the observed 443 flat-spectrum sources 
to exhibit 2 to 10% rms flux density variations in at least one epoch, shown to be dominated by 
ISS through a str ong correlation w ith Galactic latitudes as well as line-of-sight Galactic electron 



column densities (jLovell et al.ll2008l ) 



One of the most surprising results of the MASIV survey was the discove ry of a drop in IS S 
variability amplitudes and the fraction of scintillating sources at redshifts z > 2 (jLovell et al.ll2008l ) . 
Follow-up observations to the MASIV Survey at 4.9 and 8.4 GHz provided confirmation of this 
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effect at both frequencies (|Koav et al.l l2011al . hereafter Paper I). Although only a sub-sample 
consisting of 128 sources from the original survey was observed, the sources were monitored over an 
11 day period as op posed to the three or four day epochs in the MASIV Survey. We also note that 



11 aay period as op posed to tne tnree or tour aay epocns m tne iviAbiv burvey. we also note tnat 
Ofek &; Fraill ( 20 111 ) report no significant correlation between ISS amplitude at 1.4 GHz and source 



redshift, based only on 9 sources for which redshift data were available, in their investigation of the 
1.4 GHz variability of sources in the NRAO VLA Sky Survey (NVS S) and the Fain t Imag es of the 
Radio Sky at Twenty centimeters (FIRST) catalogs. Interestingly, iRichards et al.l (|201ll ) found a 



3<7 significance drop in the variability amplitudes of sources at z > 1 in their 15 GHz observations 
of about 1500 Fermi-candidate blazars over a duration of 2 years. While these longer timescale 
variations at 15 GHz are mainly intrinsic to the sources, it is conceivable that if these high-redshift 
sources are less compact (in angular and linear scales) than their low-redshift counterparts, they 
will exhibit lower levels of both ISS and intrinsic variability. 

Determining the origin of this redshift dependence of AGN ISS has important cosmological con- 
sequences, potentially allowing the ISS of AGNs to be used as a cosmological probe with the highest 
angular resolution possible. To understand this redshift dependence, we envision a simple model in 
which the flat-spectrum AGN consists of an ultra-compact, scintillating core component (estimated 
to be < 150//as from the MASIV observations) surrounded by more extended, milliarcsecond jet 
components that do not scintillate. The suppression of ISS at redshifts z > 2 can therefore be 
attributed either to an increase in the apparent angular diameters of the core components, or a de- 
crease in the compact fraction of the sources, i.e. the emission from the ultra-compact scintillating 
component of the source becomes less dominant relative to that of the extended non-scintillating 
components. Possible interpretations include: 



1. Decrease in observed brightness temperature due to cosmological expansion (IRickett et al 
20071 ) — In a sample of brig htness temperatu re-limited sources in the emission frame, as is often 



assumed for radio-loud AGNs (jReadheadlll994h . the mean observed brightness temperature decreases 
with redshift as a result of the expansion of the Universe. In a flux-limited sample, this results in 
an increase in the apparent angular diameters of sources with redshift. If this effect dominates, this 
provi des an angular size-r edshift relation for extragalactic radio sources that has long been sought 
after (|Gurvits et al.lll999h . 



2. Evolution of AGN morphologies — The drop in ISS can also be explained by a prevalence 
of sources with lower Doppler boosting factors at high redshift, which would result in either lower 
source compact fractions or larger angular diameters of the core components. However, little is 
currently known about the evolution of AGN core-jet morphologies with redshift, critical in studies 
of feedback processes in galaxy formation. It is therefore a target of very high resolution observations 
such as RadioAstron's early science programs. 

3. Source selection effects — It is well known that a flux-limited sample of fl at-spectrum AGNs 
is sev erely biased with redshift due to source orientation and relativistic beaming (ILister Marscher 
19971 ). Variations in the distribution of intrinsic source luminosities and jet Doppler boosting factors 
with redshift can lead to a redshift dependence of the apparent mean angular sizes or source compact 
fractions. Furthermore, a survey at a fixed frequency observes the sources at increasing rest frame 
emission frequencies with increasing redshift, and thereby observes the optically thick sources at 
increasing optical depths with increasing redshift, as well as different portions of the spectrum of 
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emission. 



4. S catter broadening i n the ionized intergalactic medium (IGM) and the ISM of intervening 
galaxies (IRickett et al.l 120071 ) — Cosmo-logica l models demonstrate that supernovae-driven galactic 



outflows can inject turbulence into the IGM (jEvoli &: Ferrarall201ll ). Such turbulence in the ionized 



IGM can cause angular broadening due to multipath propagation of the scattered waves. If indeed 
the redshift dependence of ISS is a result of angular broadening in the IGM, it would present a 
first direct detection of scattering in the ionized IGM, and would open up a ne w method of probing 
the p hysics of the IGM where 90% of the baryons in the Universe reside (jFukugita &; Peebles 
2004). As scattering is sensitive to the ionized components, it will complement Lyman-a studies 



of the neutral component. It may even provide an alternative means of detecting the Warm-Hot 



(Cen & Ostriker 
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it to detect ( 


Bregman 


2007 


). This scatter broadening 



may even occur in the ionized ISM of intervening galaxies, which will provide information on their 
turbulent properties, although this effect is unlikely to dominate in the majority of the high redshift 
MASIV sources. 

5. Gravitational lensing by foreground sources — The combined data from the Cosmic Lens 
All-Sky Survey (CLASS) and the Jodrell Bank VLA Astrometric Survey (J VAS) revealed only 22 
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aa. 



gravitational lens systems out of a sample of 16,503 flat-spectrum sources (jMvers et al. 
They were however, searching mainly for multiply imaged sources with arc-second resolution (with 
follow-up observations using VLB A and MERLIN for confirmation), thus would not have detected 
any low-level magnification in the sources caused by weak lensing. However, if weak lensing broadens 
the source images by 10 to 100 /xas, it could supress the ISS of sources at high redshift. Such an 
explanation would have implications for the distribution of matter (both dark and baryonic) in the 
low redshift Universe. 

Three lines of investigation are actively being pursued by the MASIV collaboration and others 
to better understand this ISS redshift dependence. The first involves a thorough examination 
of the optical properties of the MASIV sources to detect possible biases due to the presence of 
sources drawn from different AGN populations (Pursimo et al., submitted). This includes new 
observations to obtain more reliable redshift estimates and optical identifications (IDs) for the 
sources to complement archival data. New spectroscopic redshifts and IDs for many of the MASIV 
sources in which such data were not available, have also been obtained. 

The second line of investigation makes use of VLBI data to exam ine the mas-scale structures 
of the sources to determine their effects on ISS (jOiha et al.l l2004al ). While it has been found 
that scintillating s ources are more core dominated than the non-scintillating sources at mas-scales 
(lOjha et al.ll2004bl ). how their mas compact fractions scale with redshift is still being investigated. 
Multi-frequency VL BI observations to study possible angular broadening in a subsample containing 

2006) found no evidence of scatter broadening in the IGM at the 



49 MASIV sources (lOiha et a 



resolution probed (jLazio et al 



20081 ) . providing an upper limit of 500//as at 1.0 GHz. 



The third key to understanding the redshift dependence of ISS, which is the focus of this present 
paper, is to examine how this redshift dependence scales with observing frequency. The motivation 
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is that scatter broadening scales roughly with v~ 2 , whereas intrinsic source size effects scale with 
v~ l for a synchrotron self-absorbed source. On the other hand, the effects of cosmological expansion 
and gravitational lensing are achromatic. Therefore, while we gener ally expect the amplitude of 



weak ISS to decrease when we go to higher observing frequencies (see iNarayanl ((1992) for a review 
of the different regimes of ISS), we will observe either a similar or weaker redshift scaling of ISS 
amplitudes depending on which interpretation is correct. At the very least, it will allow us to rule 
out some of the interpretations listed above, or place strong constraints on them. 

The dual-frequency MASIV follow-up observations of ISS provide us with the data most suited 
for such an investigation. Indeed, the redshift dependence was found to be marginally weaker at 8.4 
GHz as compared to 4.9 GHz, which can be interpreted as due to weaker IGM scatter broadening 
of the z > 2 sources at the higher frequency (see Paper I). However, the presence of subtle selection 
effects such as those mentioned above, along with the complexity of the ionized ISM and the sources 
themselves, complicates the interpretation. 

In this follow-up paper, we present a comprehensive analysis of the data, taking these selection 
effects into consideration. We also make use of more accurate models and Monte-Carlo simulations 
to interpret the data. We provide a brief summary of the observations and characterization of source 
variability in Section [2j Section [3] then delves into the analysis and interpretation of the results. 
Our conclusions are summarized in Section SJ 



2. Observations and Source Variability Characterization 

We observed a sub-sample of 140 sources drawn from the original MASIV Survey over a duration 
of 11 days using the VLA. 70 of these sources have redshifts of z < 2 (we refer to them as the 
low-redshift sources) while another 70 have redshifts of z > 2 (the high-redshift sources). The 
telescope was divided into 2 subarrays, one observing at 4.9 GHz and the other observing at 8.4 
GHz simultaneously. Each source was observed for 1 minute at ~ 2 hour intervals. The visibilities 
were then coherently averaged over each 1 minute scan to produce the time-series data (lightcurves) 
for each source at both frequencies. 

We used the structure function, D (, s {t), as a standard characterization of the source variability 
amplitudes at various timelags r, given by: 

D obs (r) = ([S(t + T)-S(t)} 2 ), (1) 

where S(t) is the flux density of the source at time t, normalized by its mean flux density calculated 
from the entire 11-day period. The angular brackets indicate averaging over time, t. We then fit 
the following model to D i, s {t) for each source: 

7- 

-^modv ) — D sa i ■ + Z? no j se , (2) 

assuming that ISS approaches a stationary stochastic process when observed over a duration much 
longer than its characteristic timescale. D sa t is the value at which D m od(j ) saturates, and T c h ar is the 
characteristic timescale at which the structure function (SF) reaches half of its value at saturation. 
Any variability caused by instrumental and systematic errors were assumed to contribute a white 
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additive noise, Dnoisei across all time-lags. We estimated D n oise &s a quadratic sum of the flux 
independent errors i.e. system noise and confusion, and the flux dependent calibration errors, 
which we then subtracted from D moc i(T) across all time-lags to obtain the 'true' variability, D(r). 
12 sources were removed from our sample due to large errors that were not well quantified by our 
estimation of D no i se . We then used D(t) at a time-lag of 4 days, D(4d), for all subsequent analyses 
as a standard characterization of source variability amplitudes. In this paper, we denote D(4d) at 
4.9 GHz and 8.4 GHz as Z?4.9(4d) and Z?8.4(4d) respectively. Significant correlation between D(t) 
and line-of-sight Ha intensities at timelags of 1 to 7 days (see Paper I) confirm that ISS dominates 
D (t) at these timescales. 

We refer the interested reader to Paper I for an extensive elucidation of the techniques and 
analyses summarized here. The variability amplitudes of the 128 sources used in our sample are 
listed in Table [T] along with other source properties used in our analyses. 



3. Analysis and Interpretation 

3.1. Selection Effects 

The results of the MASIV Survey and the follow-up observations show that the presence of 
ISS is stron gly correlated wi th source spectral indices, mean flux densities and line-of-sight Ha 



intensities (lLovell et al.ll2008l ). We therefore examine here whether these parameters are similarly 
distributed across the low and high-redshift samples, as well as other selection effects that may bias 
the interpretation of the results. 

The sources in the original MASIV sample were selected to be flat-spectrum, core-dominated 
AGNs, based on the spectral index criterion of af'| > —0.3 (where S oc v a ). The cutoff was set at a 
higher than usual value of —0.3 (where —0.4 or —0.5 is normally used), to avoid the tail end of th e 



distribution of the classical steep-spectrum sources peaked at a ~ —0.7 (IScheuer fc Williamslll968l ). 
considering that the a\'\ values were estimated from non-coeval mean flux densities at different 
frequencies. However, we found 30 sources with spectral indices of a|'| < —0.3 in the present 
sample, of which 15 have a|;| < —0.4, calculated from the coeval mean flux densities in our 11 day 
observations. 

We also found that the redshift dependence of ISS in our sample can be at least partially linked 
to the steepening of the mean values of afg with redshift (see Figured]). The Kolmogorov-Smirnov 
(K-S) test rejects the null hypothesis, that af g of the low and high-redshift samples are drawn from 
the same parent population, with a 7.4 x 10 -3 probability that this result was obtained by chance 
(here and in all subsequent analyses, we consider probabilities below the standard value of 0.05 to be 
statistically significant, while probabilities above this value are considered not significant). Of the 
15 a|[g < —0.4 sources, which are known to scintillate less than the af;g > —0.4 sources (Paper I), 
13 are found at z > 2. Furthermore, eight of the 11 a|J > 0.4 sources are found at z < 2, and there 
are indications that their scintillation amplitudes may be larger than that of the —0.4 < a|g < 0.4 
sources (Paper I). This z-af'g correlation in itself is of great interest and we defer a full discussion 
of its significance to Section 13.41 
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Since a|;| is based on the coeval flux densities of the follow-up observations, we have used —0.4 
as a lower cutoff for the selection of sources. This allows us to retain slightly more sources for better 
statistical representation. Also, an examination of the 15 sources with —0.4 < afg < —0.3 finds 
that they are distributed roughly equally, with 8 at z < 2 and 7 at z > 2. We therefore remove only 
11 af® > 0.4 and 15 af g < —0.4 (a total of 26) sources from our sample, after which the K-S test 
shows that the distribution of afg in the high and low-redshift samples no longer differ significantly. 
In any case, we note that all subsequent analyses in this and the following sections were performed 
on both the —0.4 < afg < 0.4 and —0.3 < afg < 0.4 samples, for which we found no differences in 
the conclusions. From here onwards, we present only the results for the —0.4 < afg < 0.4 sample 
of 102 sources, comprising 46 sources at high redshift and 56 sources at low redshift. 

Considering only the —0.4 < af'| < 0.4 sources, we found that the drop in ISS amplitudes at 
high redshift remains statistically significant. The K-S test confirms that the variability amplitudes 
of the high-redshift sources, characterized by D^g(4d) and Z?8.4(4d), are significantly smaller than 
that of the low-redshift sources, with a probability of 5.1 x 10 -5 that this occured by chance. This 
is clearly seen in the histograms of Figure [2j 

We note that in the selection of sources for the MASIV follow-up observations, the z < 2 sample 
was biased towards the variable sources, while all sources with known redshifts were selected for the 
z > 2 sample. Although one could argue that the fraction of scintillating sou rces at low redshif t 



is higher than that at high redshift anyway, as found in the MASIV survey (jLovell et al.l 12008). 
this introduces another possible selection effect into the source sample. However, the significant 
decrease in ISS amplitudes at high redshift of the —0.4 < afg < 0.4 sources cannot be attributed 
solely to this selection effect, as the significance is retained even when only the most variable 
sources are selected. K-S tests show that -D4.g(4d) and -Ds.4(4d) in the z < 2 sources are still 
significantly larger than their counterparts in the z > 2 sources when considering only the 72 
sources where D(4d) > 2 x D no i se (equivalent to > 3a variability amplitudes) at both frequencies, 
with probabilities of 1.3 x 10 -2 and 2.2 x 10~ 2 that they occured by chance. These 72 sources were 
drawn from only the —0.4 < afg < 0.4 sample. 

The sources in our sample were also carefully selected so that their mean flux densities and line- 
of-sight Ha intensities would be evenly distributed in the high and low-redshift samples, but a large 
fraction of these sources are variable, so their mean flux densities may have changed. As a check, 
we performed the K-S test which found no statistically significant differences in the distribution of 
mean flux densities and line-of-sight Ha intensities in the low and high-redshift samples. This is 
true before and after the removal of the 26 afg > 0.4 and afg < —0.4 sources, as well as for the 72 
sources with > 3er variability. 

Another possible source of selection bias is the presence of sources with different optical IDs in 
the sample. 45 out of the 46 high-redshift sources are identified as flat-spectrum radio-loud quasars 
(FSRQs), while the other source has yet to be identified, most likely an FSRQ. In the low-redshift 
sample however, 42 sources are identified as FSRQ s, 12 are BL Lac obj ects, and 2 are Seyfert 
1 galaxies. In standard AGN unification schemes ( Urry &: Padovanilll995l ). FSRQs and Seyfert 1 
galaxies are classified as Type 1 AGNs, with both having broad emission lines. They differ only 
in their radio and optical luminosities, the significance of which is historical rather than physical, 
so we group them together in our analysis. On the other hand, BL Lacs are classified as Type 
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AGNs due to their weak or lack of emission lines. BL Lacs have been observed to be more variable 



than FSRQs, intrinsically (IRichards et al.ll201ll ) as well as in terms of their ISS (Pursimo et al., 
submitted). We therefore carry out our analysis with and without the inclusion of the BL Lacs to 
determine if they affect the interpretation of the data. 

While biases caused by the aforementioned parameters can be mitigated through the careful 
selection of sources, there are other selection effects that are unavoidable in a survey such as this. 
These selection effects can increase or decrease ISS with redshift, biasing the result either way. 
For example, the sources are observed at increasing rest-frame emission frequency with increasing 
redshift. For an optically thick synchrotron self-absorbed source, the source size decreases with 
increasing rest-frame frequency. 

Furthermore, a flux-limited survey will always be affected by the Malmquist bias arising from 
the scaling of source spectral luminosity, L^, with redshift. In fact, this perceived redshift depen- 
dence of ISS may even be interpreted as a luminosity dependence. It is possible that the higher- 
luminosity sources are intrinsically larger, and may therefore scintillate less. The plot of the 4.9 
GHz spectral luminosities (calculated using Ho = 70 kms _1 Mpc _1 , Qm = 0.27, Qa = 0.73 and as- 
suming isotropic emission) against source redshifts (Figure [3]) shows that at any particular redshift, 
there is a luminosity dependence of ISS because the stronger sources (S4.9 > 0.3 Jy) scintillate less 
than the weaker sources (S4.9 < 0.3 Jy). There is a visible redshift cutoff of Z?4.9(4d) at z > 2 for 
both the weak and strong sources. There are also possible luminosity cutoffs at L4.9 ~ 10 28 WHz _1 
and L4.9 ~ 10 27 WHz _1 for the strong and weak sources respectively. However, it is difficult in 
reality to determine if the ISS amplitudes of AGNs (in our present sample at least) have a redshift 
cutoff, luminosity cutoff or both, since S u , z and L u are inter-dependent. This problem is further 
compounded for relativistically beamed sources, where Doppler boosting complicates the definition 
of the intrinsic source luminosities. Intrinsically more luminous sources can be detected at larger 
angles of orientation to the line of sight, and may thus have lower Doppler boosting factors which 
decrease their compact fractions with redshift. On the other hand, if the high-redshift sources 
are dominated by highly Doppler-boosted sources with intrinsic luminosities comparable to their 
low-redshift counterparts, the compact fractions of the high-redshift sources may be larger. 



3.2. Structure Function Ratios 

We present a method to alleviate the effects of the selection biases discussed in Section 13. 1\ 
in particular the redshift dependence of source luminosities, Doppler boosting factors and rest- 
frame emission frequencies. Instead of comparing the redshift dependences of the mean values 
of -D4.g(4d) and -Dg.4(4d) separately, as was done in Paper I, it is more instructive to obtain the 
ratio of Z?8.4(4d) to Z?4.g(4d) for each individual source, then compare the mean and distribution 
of the ratios at high and low redshift. This normalizes each source by itself, thereby reducing 
the dependence of the results on these source parameters if they are frequency independent or 
have known frequency dependences. The rest-frame emission frequencies also become irrelevant, 
since we are only interested in the ratio of the two observing frequencies. Parameters involving 
the properties of the ISM are also factored out, again assuming all these source parameters to be 
frequency independent. We provide a mathematical justification for these claims in Section 13.2.11 
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and present our results in Section [3.2.21 



3.2.1. Theoretical Basis 



To obtain theoretical estimates of the structure function ratios for comparison with the obser- 
vational data, we make use of standard ISS models in which the scattering region is approximated 
as a thin, phase-changing screen with an isotropic Kolmogorov spectrum. The quantity of interest 
is the spatial coherence, T±{r;v), of the flux measured a t two l ocations separated by a distance r 
on the Earth at a frequency of v ([Goodman fe Naravanl (|1989l ) provide the detailed formalisms). 
The model assumes that the phase structure function at the scattering screen is isotropic at the 
length-scales of interest and does not vary at the timescales (r) of interest, so that the spatial 
coherence can be simply related to the temporal coherence as measured by a single telescope by 
equating r = v s t, where v s is the transverse velocity of the scattering screen relative to the Earth. 
We then compute £>(4d) = 2[r 4 (0; u) - T 4 (r = v s ■ 4d; v)}. 

Anal ytical solutions for Tj{r\ u) in th e very weak and very strong scintillation regimes are 



given m 



Coles et al. I (j 19871 ) and iNaravanl (|1992l ). They use the modulation index, m, defined 



as the rms variations as a percentage of the mean flux density of the source, to quantify the 
variability amplitude of the source. Assuming that the structure functions saturate within 4 days, 
D(4d) ~ 2m 2 . In weak scintillation (y S> vt, where vt is the transition frequency bet ween weak an d 
strong ISS), the modulation index of a point source is given in the following form by I Walker! ((1998): 



m 



p,w 



U t \ 17/12 
V J 



(3) 



For strong refractive scintillation [y <C ff), the point source modulation index is ([ Walker! 1 19981 ): 

17/30 



p,r 



(4) 



The modulation index of an extended source is then further suppressed relative to that of a point 
source by a factor tha t is dependen t on the apparent angular size of the source, 9, as it appears to 
the scattering screen ()Walkerlll998l ): 



m 



'ext 



7/6 



(5) 



where m p is equivalent to m p ^ w or m Ptr . 9 ex t is the angular size above which the source can no longer 
be considered a point source. For weak ISS, 6 ext is equivalent to the angular size of the first Fresnel 



zone at the scattering screen, given by Of = y c/ (27ti/Dism) 



oc v 



-0.5 



where Dism is the distance 



from the Earth to the scattering screen and c is the speed of light. For strong refractive ISS, 9 ex t is 
the refractive scale of the density inhomogeneities at the scattering screen, given by 9 r oc v 



-2.2 



The ratio of Ds.4(4d) to Z?4.g(4d), which we denote as Rd, can be calculated using the asymp- 
totic limits of Equations [3] to [5] to obtain: 



Rd 



m 



8.4 



0.217 



m 



1.9 



Of. 



4.9 



7/3 



V#8.4 



7/3 



(6) 
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for an extended source in the weak ISS regime, and: 



7/3 //, N 7/3 
7 r,8.4 \ / ^4.9 



fl D « 1.842 p p , (7) 

for an extended source in the strong refractive ISS regime. Any anisotropic properties of the ISM 
(i.e. the strength of turbulence and distance to the scattering screen), encapsulated in the terms u t , 
Of and 9 r , either cancel out or have known frequency dependences. Rd is therefore sensitive only 
to the ratio #4.9/6*8.4- 

The apparent source size, whose frequency scaling is dependent upon whether it is dominated 
by intrinsic effects or scatter broadening, can be modeled as: 



ei c + 9 2 scat <xv-P, (8) 

where 9 src is the intrinsic source angular size, and 6 SC at represents additional contributions due to 
scatter broadening in the ISM or the IGM. If intrinsic source size effects dominate, 9 ~ 9 src , and 
any source dependent parameters that 9 src is a function of, such as the luminosity, compact fraction 
and Doppler boosting factor, cancel out in Equations [6] and [7| assuming that they are frequency 
independent. Rd is therefore transparent to any redshift dependences of these parameters. The 
ratio of the emission-frame frequencies is a constant across all redshifts for a fixed pair of observing 
frequencies, so Rd is also insensitive to source properties that vary with increasing emission-frame 
frequencies. Typically, j3 ~ 1 for a synchrotron self-absorbed source, while j3 ~ if the source size 
is frequency independent. On the other hand, /3 ~ 2.2 for a scattering screen with Kolmogorov 
turbulence, if the source size is dominated by scatter broadening so that 9 ~ 9 scat- 

We calculated Rd for different values of j3 in the asymptotically weak and strong refractive 
ISS regimes. In the case of a point source (9 < 9 ex t), Rd ~ 0.2 and Rd ~ 1.8 in the weak ISS and 
strong refractive ISS regimes respectively, and is independent of j3. At our observing frequencies 
and for sight-lines through mid-Galactic latitudes however, 9 > 9 ex t for AGN in general. In this 
case, Rd < 0.4 when intrinsic source size effects dominate (/3 < 1) and Rd ~ 1-8 when scatter 
broadening dominates (f3 = 2.2), for both weak and strong ISS. 

In the intermediate scattering regimes typical of our observations, where there are no analytica l 



solutions, we make use of the fitting functions for T^r; v) provided bv lGoodman fc Naravanl (|2006l ) 
based on numerical computations, which allow us to calculate Rd when 4.9 GHz < v% < 8.4 GHz. 
Figure U] demonstrates how Rd varies with v% for f3 = 0, 1 and 2.2. 

While the analytical approach provides a better understanding of the physics involved in the 
analysis of the SF ratios, it assumes asymptotically weak or strong refractive ISS of the sources. 
Furthermore, it assumes that the characteristic timescales of ISS are less than four days so that 
Z)(4d) is well approximated by 2m 2 . The ISS timescales can vary with observing frequency, which 
affects Rd if the SFs have yet to saturate at one or both frequencies. We know from Paper I 
that ~ 20% of the sourc es have ISS timescales of mo re than four days on at least one frequency. 



On the other hand, the iGoodman fc Naravanl (|2006l ) fitting function makes no such assumption 



about the ISS timescales and simply calculates T^Vs ■ 4d; u). However, when the ISS timescales are 
longer than four days, we note that Rd estimated from the fitting function becomes sensitive to the 
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various scattering screen and source parameters, which will be true of the observations as well. In 
interpreting our data, we mainly used the fitting functions for comparisons with the observations. 



3.2.2. Observational Results 

The histogram of Rd as obtained from our observations exhibits a well defined peak in the 
0.25 < Rd < 0.50 bin (Figure [5]), consistent with the weak ISS of sources dominated by intrinsic 
source size effects so that /3 < 1 (see Figure SJ). As Rd is inaccurate for sources whose D(4d) is 
close to or lower than the noise floor at one or both frequencies, the histogram includes only the 72 
sources with > 3<r variability amplitudes at both frequencies, selected from the —0.4 < a|;g < 0.4 
sample, with 48 sources at low redshift and 24 at high redshift. Of the 48 low redshift sources, 11 
are Type AGNs (BL Lacs) and 37 are Type 1 AGNs (FSRQs and Seyfert 1 galaxies). 

However, the tail towards larger values of Rd indicates the presence of at least another im- 
portant effect. Three effects potentially increase Rd so that it becomes comparable to or greater 
than unity. One is if vt > 5 GHz, so that the sources are scintillating in the strong ISS regime or 
at the transition between weak and strong ISS at one or both frequencies. The second possibility 
is that (5 > 1 due to scatter broadening at a second, more distant screen before the waves arrive 
at the scintillation-inducing screen. This second scattering screen can be Galactic or extragalactic. 
The third possibility is that the SFs have yet to saturate within 4 days, so that the assumption of 
Z)(4d) ~ 2m 2 no longer holds. Since r c h a r oc 6, and assuming 9 oc f , the scintillation timescales 
are shorter at 8.4 GHz. This causes the 8.4 GHz SFs to rise faster and saturate earlier in comparison 
to that at 4.9 GHz, thereby increasing Rd- 

We determined that the sources with > 1 are not significantly scatter broadened, based on 
an examination of the Rd values in the weak (£4.9 < 0.3 Jy) and strong (S4.9 > 0.3 Jy) sample of 
sources. We can assume that the weak and strong sources have similar mean intrinsic brightness 
temperatures, so that a ~ 0.1 Jy source tends to have a smaller angular diameter than a ~ 1.0 Jy 
source. This is not an unreasonable assumption, as it explains why the weak sources have been 



observed to scintillate more than the stronger sources (jLovell et al.l 12008). Additionally, the lower 
ISS amplitudes in the strong sample of sources cannot be attributed to the presence of stronger 
mas-scale jet components, as confirmed by VLBI observations that found no significant diffe rence in 



the mas-scale morphologies of the strong and weak flux density sources (jOiha et al.ll2004bl ). From 
Equation [81 we see that the effects of scatter broadening will be more dominant in sources with 
smaller intrinsic angular sizes. Therefore, if the sources in our sample are scatter broadened, we 
should observe higher values of Rd in the weaker sources than in the stronger sources. K-S tests do 
not show that the weak sources have Rd values significantly higher than that of the strong sources. 
This is true when all 72 sources are considered and when only the Type 1 AGNs are considered. 
In fact, the converse is true. Figure [6] shows scatter plots of -Dg.4(4d) against D^ g(4d) for the 
72 sources with D(4d) > 3a, classifying the sources into weak and strong samples. The dotted 
lines have slopes of 0.4 and 1.8, representing the possible values of Rd in the asymptotic weak and 
strong ISS regimes. The solid lines represent linear least-square fits for the two source categories. 
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In obtaining the fits, each source is weighted by a factor: 

w = (°D,4.9 + °D,&.i)~ ' > ( 9 ) 

where an, 4.9 and (?d,%a are the normalized errors in -D4.9 (4d) and L>8.4(4d) respectively. This means 
that sources that have smaller errors in D(4d) have larger weights in the fitting process. The dashed 
lines are the 99% confidence bounds for those fits. It can be seen that the strong sources have a 
larger fitted Rd of 1.06 ± 0.34 as compared to a fitted Rd of 0.44 ± 0.18 for the weak sources at 
Vo confidence. 



We found Rd to be influenced by the strength of ISS as indicated by the line-of-sight Ha 
intensit y to each source, ob tained from the Wisconsin H-Alpha Mapper (WHAM) Northern Sky 



Survey (jHaffner et al.ll2003l ). The line-of-sight Ha intensity is in units of Rayleighs (R), and we 
denote it as I a . The K-S test confirms that Rd in the I a > 5 R sample is significantly larger than 
that of the I a < 5.0 R sample, with a 1.3 x 10~ 2 probability that this occured by chance when all 
sources are considered, and a probability of 4.0 x 10 -2 when only the Type 1 AGNs are considered. 
The scatter plots and fits of Rd in Figure [7] also show that sources with I a > 5 R tend to have 
larger values of Rd- The fitted Rd is found to be 0.99 ± 0.59 for I a > 5 R as compared to 0.35 
± 0.08 for I a < 5.0 R at 99% confidence. The sight-lines with larger Ha intensities have higher 
electron column densities and are at lower Galactic latitudes where the sources are seen through 
thicker regions of the Galaxy. Therefore, we expect the transition frequencies to be higher through 
these sight-lines. This demonstrates that sources with Rd > 1 are scintillating in the strong ISS 
regime or at the transition between weak and strong ISS. Additionally, the five most variable sources 
with Rd > 1 all have I a > 5 R, consistent with scintillation amplitudes being the highest at the 
transition frequency between weak and strong ISS. 

In comparing the values of Rd in the z < 2 and z > 2 subsamples, we rule out scatter 
broadening in the IGM as the origin of the redshift dependence of ISS, in our present sample at 
least. Using the K-S test, we find that Rd in the high-redshift sample is not significantly higher 
than Rd in the low-redshift sample considering all sources and only the Type 1 AGNs. Figure [8] 
shows only a marginal increase in the fitted Rd of sources in the z > 2 sample as compared to the 
z < 2 sample, with 0.83 ± 0.35 for z > 2 and 0.41 ± 0.16 for z < 2 (99% confidence intervals). 
The fitted Rd values at low and high redshift are calculated separately for the weak and strong 
sample of sources, summarized in Table El The redshift dependence of Rd in the weak sources is 
not significantly steeper than that of the strong sources, providing more evidence against scatter 
broadening in the IGM or that it has any redshift dependence at the resolution of our observations. 

In fact, the marginal increase in Rd with redshift is most plausibly attributed to intrinsic 
source size effects, since the stronger sources that have larger intrinsic angular sizes also have larger 
fitted Rd- We explain this increase in Rd with redshift in Section \3.3\ and discuss the implications 
of these results on IGM scattering properties in Section 13.51 
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3.3. Decrease in Observed Brightness Temperature Due to Cosmological Expansion 



We propose that the suppression of ISS at z > 2, considering only the —0.4 < af;g < 0.4 
sources, can be sufficiently explained by the decrease in observed brightness temperature of a flux- 
limited sample of sources due to cosmological expansion. The angular size of a source, 6 src is related 
to its observed brightness temperature, Tb t0 b s through the following well-known expression: 



r 2 S" 

(10) 



b,obs 

where v is the observing frequency, S v is the observed flux density, c is the speed of light and k is 
the Boltzmann constant. On cosmological scales, the observed brightness temperature is a factor of 
(1 + 2:) lower than the brightness temperature in the frame of emission, Tb em , due to the expansion 
of the Universe. em is in turn a function of the intrinsic brightness temperature of the source, 
Tb,int, Doppler-boosted by a factor 5 due to the effects of relativistic beaming in AGN jets. Since 
only the most compact components of the source core scintillate in the presence of turbulence in the 
ISM, and it is the sizes of these compact regions that we are inferring from the source variability, 
we multiply S u by the compact fraction, f c of the source. Equation [TO] then becomes: 



{\ + z)<?fJ3 v 

This is the angular diameter of the source as it appears to the scattering screen in the ISM, assuming 
no additional increase in apparent size due to extrinsic propagation effects. Therefore, 9 src cx 
(1 + z) 5 and we can expect the ISS amplitudes to decrease with redshift, contingent upon the 
following assumptions: 

1. The mean flux densities of the sources do not vary with redshift, i.e. the sample is flux-limited. 

2. The intrinsic brig htness temperatures, Th^nt, have a cutoff, either at the inverse Compton 



limit at ~ 10 12 K (IKellermann fc Paulinv-Tothlll969l ). or at the energy equipartition limit at 
~ 10 11 K dReadheadlll994i : lLahteenmaki et alJll999h . 



3. The mean Doppler boosting factors and compact fractions of the sources remain constant and 
do not evolve with redshift. 

4. Any possible effect of decreasing angular sizes of the optically thick cores with increasing 
rest-frame emission frequencies is ignored. 

Based on this model, we performed numerical computations using the fitting function in 
Goodman fc Narayan (j20od ). plotting L>(4d) of the weak and strong sources against redshift for 



various values of the Doppler boosting factor in Figure [H We used the fiducial values shown in 
Table [3] as the model parameters. We assumed that the phase fluctuations at the scattering screen 
obey a power law with a Kolmogorov spectrum, and that the sources have a Gaussian intensity 
profile. The mean D(4d) obtained from the observations are shown for two redshift bins, separating 
the weak (54.9 < 0.3 Jy) and strong (S4.9 > 0.3 Jy) sources. 
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The observed redshift dependence of the mean values of Z?(4d) in the 102 -0.4 < a|;| < 0.4 
sources appears to be consistent with the model and its assumptions. Even with the possible bias 
towards more variable sources at z < 2, the decrease in mean ISS amplitudes is no greater than 
that expected from this model. As further confirmation, this agreement holds, within the la error 
bars, even when only the 72 > 3a variable sources were used in obtaining the mean values of Z?(4d). 
Looking at Figure [9j we see that the model successfully explains the weaker redshift dependence 
of Z?8.4(4d) as compared to that of Z?4.g(4d), without having to invoke scatter broadening in the 
IGM. Furthermore, it explains why the ISS amplitudes of the strong sources have a weaker redshift 
dependence than that of the weak sources, as also reported in Pursimo et al. (submitted) for the 
larger MASIV sample. 

The observations and assumptions of the model, particularly that of constant mean Doppler 
boosting factors with increasing redshift, are also con sistent with other statistic al studies of Doppler- 
boosted AG N sources. Monte C arlo simulations bv lLister &: Marscherl (j 19971 ) and recent observa- 
tional data (jHovatta et al .1120091 ) suggest that the mean Doppler factors for a flux-limited sample of 
sources remains relatively constant at z > 0.6. The highly beamed sources with large Doppler fac- 
tors are very rare, since their jets need to be aligned very close to the line of sight to be detectable. 
Considering that these large S sources would also appear very luminous, one would expect to detect 
more of these sources with increasing redshift where the volume is also larger (thus increasing the 
likelihood of detecting these rare sources). However, this is offset by the large range in intrinsic 
luminosities of the sources, so that sources with large intrinsic luminosities and low Doppler factors 
(due to large angles of orientation) will also be included in the high-redshift sample. While the 
range of source Dopple r factors increases with re dshif t in a flux-limited samp le, the mean remains 
the same. According to lLister Sz Marscherl ( 19971 ) and lArshakian et al.l (201CJ), at redshifts z < 0.6, 



the mean Doppler boosting factor is lower due to the scarcity of sources with large Doppler factors 
(5 > 20) within the small volume at such redshifts. This may explain why the most variable sources 
are not found below z < 0.6 at both frequencies, as can be seen i n the scatter plots of Figure M 
This effect is seen in the original MASIV data as well (Figure 13 in lLovell et al.l ()2008l )). 



The marginal increase in Rd with redshift as seen in Figure[8]and Table[2]can also be sufficiently 
explained with the same model of decreasing Th,nhs with redshift. Figure [101 shows model values of 
Rd (blue and red curves), calculated using the lGoodman Naravanl (|2006l ) fitting formula and the 
same input parameters in Table O The binned plots depict the fitted Rd at low and high redshifts, 
with the error bars given by the 68% confidence bounds (roughly equivalent to la errors). The 
model calculations show that Rd can indeed increase with redshift without including any scatter 
broadening effects. This redshift dependence of Rd arises due to the increase in source angular size 
with redshift, which in turn increases the timescale of the scintillations for a fixed scattering screen 
velocity. As explained in Section 13.21 the timescales can increase sufficiently so that the SFs do not 
saturate within four days, leading to higher values of Rd- This also explains why the sources with 
high Rd also tend to be strong sources with larger angular sizes rather than weak sources, as seen 
for the fitted Rd from the observations. To further illustrate this point, we include in Figure PTOl the 
model values of Rd for the case where the scattering screen velocity is reduced to 20 kms -1 for the 
strong sources, shown as black curves. In this case, the SFs take a longer time to saturate, thereby 
increasing Rd even further. While the increase in fitted Rd from low to high redshift appears larger 
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than that predicted by the model for the weak sources, this difference is < 2a. In any case, this 
difference cannot be attributed to scatter broadening since the fitted Rd of the strong sources is 
clearly larger than that of the weak sources at both low and high redshift. 



We applied the Monte Carlo method to the iGoodman <fe Naravanl (|2006l ) fitting functions to 



obtain simulated distributions of I?4.g(4d), Z?s.4(4d) and Rd of each source for further comparisons. 
Figure [TT] shows observed values of -D4.g(4d), -Ds.4(4d) and Rd of the —0.4 < all < 0.4 sources (72 
sources with > 3a variability in the case of Rd) plotted against their corresponding median simu- 
lated values. The horizontal error bars represent the median absolute deviations of the simulated 
values of Z)4.g(4d), _Ds.4(4d) and Rd for each source. 

For each source, the 4.9 GHz and 8.4 GHz mean flux densities, as well as source redshift, 
are kept constant at their observed values, since these parameters are known. We also keep the 
transition frequency and scattering screen distance of each source constant, estimated from the 
line-of-sight Ha intensity and Galactic latitude of the source (see Appendix [A] for details). We fixed 
the intrinsic brightness temperatures and scattering screen velocities at the typical values of 10 11 K 
and 50 kms -1 respectively for all sources. We then randomly generated 1000 values of the Doppler 
boosting factor and source compact fraction for each of the sources, with Gaussian distributions 
peaked at 15 and 0.5, and standard deviations of 4 and 0.1 respectively. 

The observed values of logio[-D4.g(4d)], logio[£>8.4(4d)] and logiof-Ro] show statistically sig- 
nificant correlations with their simulated counterparts. We obtained Pearsons linear correlation 
coefficients of 0.54 for both logio[D4.g(4d)] and logio[-Ds.4(4d)] respectively, with probabilities of 

4.0 x 10~ 9 and 3.9 x 10~ 9 that they were obtained by chance. We obtained a weaker correlation 
coefficient of 0.23 for log[i?o], with a probability of 0.05 that this was obtained by chance. The 
correlation is weaker for Rd due to the presence of sources with observed Rd values below 0.4, 
possibly due to errors in the estimation of Rd or intrinsic source sizes with frequency dependences 
flatter than the used in our model. 

Our simulated dataset exhibits similar trends to that of the observations. K-S tests performed 
on the simulated SFs show that D^g(4d) and Z?8.4(4d) of the z < 2 sources are larger than that 
of the z > 2 sources, with probabilities of 6.1 x 10 -5 and 9.0 x 10 -5 that these were obtained by 
chance. As in our observational data, the high-redshift Rd values are not significantly larger than 
the low-redshift Rd values, even though the six sources with the largest simulated Rd values are 
all at z > 2 (Figure fTT]) . Additionally, the K-S tests show that Rd values for both the £4. 9 > 0.3 Jy 
and I a > 5.0 R source samples are larger than Rd values in the £4. 9 < 0.3 Jy (with a probability 
of 2.8 x 10" 4 that this result was obtained by chance) and I a < 5.0 R (with a probability of 

1.1 x 10 -3 that this result was obtained by chance) source samples, broadly consistent with that of 
our observations. 

This model ofD(r) oc (1+z) 5 is consistent with that of the original 4.9 GHz MASIV dataset as 
well, for ~ 250 sources in which redshift data are available, except at z > 3 where the observed D(2d) 
is ~ 2a below the predicted curves (Pursimo et al., submitted). However, this steeper than predicted 
drop in D(2d) can be attributed to the additional effect of the z-a correlation, which cannot be 
further analyzed for the original MASIV sample due to observations at only a single frequency. The 
model and Monte Carlo simulations will hence need to be applied to the full MASIV dataset to 
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further test this model of decreasing observed brightness temperatures, pending the acquisition of 
optical spectroscopic redshifts for the remaining sources for which they are not available. The larger 
sample may also allow us to break the redshift-luminosity degeneracy, providing a stronger test of 
whether we are observing a redshift cutoff or luminosity cutoff in D(t). 



3.4. The z-ck^g Correlation: Selection Effect or Source Evolution? 



We discuss why the mean source spectral indices steepen with increasing redshift (Figure [T]), 
and why 13 out of the 15 sources with afg < —0.4 in the sample lie at z > 2. This effect partially 
accounts for the redshift dependence of ISS in our sample, as it has been established in Paper I that 
the mean Z?(4d) of the a|'g < —0.4 sources is a factor of 



6 lower than that of the a|g > 



-0.4 



sources. The afi < -0.4 sources therefore either have source sizes that are on average a factor of 
~ 2 larger or compact fractions that are a factor of ~ 2.5 lower than their af g > —0.4 counterparts. 

This z-a%\ correlation could simply be a selection effect. As mentioned in Section 13.11 the 
original MASIV Survey sources were selected to have a\\ > —0.3, but the sample could have been 
contaminated by sources whose spectral indices were inaccurately estimated due to the non-coeval, 
variable flux densities. The comoving spatial density of fl at-spectrum radio sources has been foun d 
to decrease by a factor of ~ 5 between redshifts 2 < z < 4 (IPeacock II 19851 : iDunlop k, Peacocklll990l ) , 
perceived as evidence for a peak in quasar activity at z y 2.5. A similar but slower decline was 
found for steep-spectrum sources (IDunlop &: Peacock! 1 199G ) . which may explain why there are more 



of these steeper af g sources at high redshift. Two of the high-redshift, steep-spectrum sources have 
been optically identified as narrow-line radio galaxies (Type 2 AGNs with narrow emission lines), 
and so may be classical steep-spectrum sources. After removing these two sources from the original 
sample of 128 sources, the K-S Test still finds a statistically significant difference in the spectral 
indices in the low and high-redshift source samples. The fact that the z < 2 sample is biased 
towards variable sources, while all sources with known redshifts at z > 2 were selected regardless 
of variability, could also have biased the low-redshift sample towards flatter af'g, and vice versa for 
the high-redshift sample. 

It is also possible that some of these sources have convex spectra that steepen at higher fre- 
quencies, as found in gigahertz peaked-spectrum (GPS) sources; k-correction effects then lead to a 
steep ening of spectral indices d ue to increasing rest-frame emission frequencies with increasing red- 
shift. iJarvis Rawlingsl (|2000l ) have demonstrated that a significant portion of the most luminous 
radio-selected flat-spectrum sources are GPS sources. Since these GPS sou rces are not D o ppler - 
boosted, they are less compact and therefore scintillate less. Furthermore, IChhetri et al.l (120121 ) 
recently found a steepening of source spectral indices with redshift in compact radio sources from 
the Australia Telescope 20 GHz (AT20G) Survey. This z-a correlation was discovered to be more 
prominent when a|;| is used as compared to af;®, which they attribute to spectral curvature and 
k-correction effects. 



A natural and physical explanation for this z-a|'| correlation is also conceivable. The steep- 
ening of source spectral indices with redshift has long been o bserved in radio galaxies (classical 



steep-spectrum sources) (jLaing fc Peacock 



1980 



Macklin Ill982l ). where spectral index cut-offs have 
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been successfully used to search for radio galaxies at the highest redshifts. Traditional expla- 
nations for th is correlation includ e (1) k-correction effects in sources with convex spectral en- 
ergy densities (|Gopal-Krishnalll988l ) (2) inverse Compton losses associat ed with Cosmic Micro wave 
Background (CMB) photons whose energy densities scale with (1 + z) 4 ( Krolik fe Chen 1991) and 



(3) a luminosity-spe ctral index relati on coupled wi t h the Malmquist bias (jLaing h Peacock 



Blundell et al.lll999l ). More recently, iKlamer et al.l (|2006l ) argue that this z-a correlation could be 



1980; 



due to higher ambient densities at high redshift, resulting in increased synchrotron and inverse 
Compton losses in pressure-confined radio lobes. The evidence comes from the observed properties 
of high-redshift radio galaxies, which include (1) similarities to low-redshift radio galaxies residing 
in dense clusters, (2) extreme rotation measures (RMs), and (3) knotty radio emission interpreted 
as frustrated jets in dense and clumpy environments. For radio galaxies, one would expect to ob- 
serve increasingly compact sources at high redshift. All known radio- loud AGNs at z > 4 are steep 



(Frev et al. 


2008, 


2010, 


2011 



thought to be young and 



'frustrated' radio galaxies. 



If the flat-spectrum, Doppler-boosted population of radio-loud AGNs reside in similar environ- 
ments at high redshift, it is likely that pressure from the surrounding IGM, whose densities scale 
with (1 + z) 3 , will reduce their Doppler boosting factors. This in turn reduces the compact fractions 
of the sources and reduces their scintillation amplitudes. This will also result in a steepening of 
spectral indices, as the contribution of the optically thick core components to the mean observed 
flux densities is reduced relative to that of the optically thin jet components. 

Multifrequency VLBI studies based on new observations and archival data will be needed to 
determine if this z-a|' 4 correlation and its relationship to the redshift dependence of ISS, is mainly 
due to selection effects or interesting physical phenomena related to AGN evolution. 



3.5. Constraints on IGM Scattering and Turbulence 

While our observations provide no clear detection of scatter broadening in the IGM or that 
it has any significant redshift dependence between the z < 2 and z > 2 subsamples, we can place 
strong constraints on it. The top panel of Figure [12] shows estimated app arent angular sizes of 



all sou rces at 4.9 GHz, #4.9, calculated from their D^g(Ad) values using the iGoodman k. Naravan 



([20061 ) fitting functions. For all sources in which Z?4.g(4d) > D no i se , the upper limits of #4.9 are 
calculated using f c ~ 1, while the lower limits are calculated using f c ~ 0.1. For sources in which 
Z?4.g(4d) < D no i se , upper limits to 6*4.9 cannot be obtained, while the lower limits are calculated using 
f c ~ 0.1 and setting D no i se as the upper limit to the source variability. For each source, we used u t 
and -Dism calculated in AppendixE] which were also used in the Monte Carlo simulations described 
in Section 13.31 The estimation of #4.9 makes no assumptions about the brightness temperatures 
and Doppler boosting factors of the sources. The upper limits of #4.9 are also effectively upper 
limits of the 4.9 GHz sca t and 8 src for the sight-lines to our sources (Equation [HJ) ; they are shown 
proportional to the sizes of the circles in the lower panel of Figure [12] in Galactic coordinates. 

We go one step further in constraining IGM scatter broadening by making use of the fitted Rjj 
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obtained for the > 3a variable sources. At 4.9 GHz, the upper limit to scatter broadening can be 
formulated from Equations [6] and [8] in the weak ISS regime as: 



'scat 



< 



src(max) 



2.l§{R Dlymax) ) 7 - 1 



l-0.59(i? D 



(max 



0.5 



(12) 



for 0.4 < RD(max) < 1-8, where RD(max) is the upper limit of Rd and sr c(max) is the upper limit 
to the intrinsic source angular size at 4.9 GHz. Again, we made use of the relations 9 src oc u^ 1 
and 6 scat ^ _2 ' 2 - This inequality posits that sources scintillating in the weak ISS regime will have 
Rd ~ 0.4 when completely dominated by intrinsic source size effects. As an increase in IGM scatter 
broadening increases Rd, the upper limit to 9 SC at is determined by the level of increase in Rd above 
this nominal value. The dominance of 9 sca t is also dependent on Ogre] sources with smaller intrinsic 
angular sizes are more likely to be dominated by 6 SC at than sources with larger angular sizes. 



We use the weak ISS approximation by iNarayanl (119921 ) and IWalkerl (|1998l ) here because it 
gives the most conservative upper limit to sc atter broadening in th e IGM and is not dependent on 
any other model parameters. Comparing the I Goodman &: Narayanl (120061 ) fitting function with the 



weak ISS model ONaravan 



1992 



Walkerlll998l ) in calculating Rd, we found that both provide similar 

As opposed t o the weak ISS model 



constraints when 9 = y 9 2 src + 9 2 cat is < 50/ias (see Fi gure fl~3l) 
where D(4d) ~ 2m 2 is assumed, Rd rises faster in the lGoodman fe Naravanl ( 20061 ) fitting function 
with increasing 9 when 9 is of an order ~ lOO^as, as the scintillation timescales exceed 4 days and 
the SFs do n ot saturate (the same r e ason why Rd increases marginally with redshift in Figure [10]) . 
Although the lGoodman &: Naravanl ((2006) fitting functions give stronger constraints, they are very 
dependent on the parameters of the model when 9 is large. 

Since the upper limit of Rd for the low-redshift sources is 0.57 at 99% confidence, we obtain 



'scat 



< 110/ias at 4.9 GHz, assuming that the intrinsic angular sizes of the scintillating components 
in all our sources are < 150/ias (as seen for the > 3a variable sources at low redshift in the 
upper panel of Figure PT2]) . In discussi ng the feasi b ility of usi ng the Square Kilom etre Array to 
detect intergalactic scatter broadening, lLazio et al.l (|2004l ) and iGodfrev et al.l (|201ll ) propose that 
angular resolutions better than 4 mas at 1.4 GHz and 80 mas at 0.33 GHz are required. Our most 
conservative constraints push these limits lower. A simple extrapolation gives 9 sca t < 1.7 mas at 
1.4 GHz and 9 scat < 42 mas at 0.33 GHz. 

The strongest constraints can be derived from the Rd of the weak sources, which is also no 
more than 0.57 at 99% confidence for the weak, low redshift sources. With flux densities of ~ 0.1 
Jy 9 src can be as low as ~ 10/uas. For example, the source PKS 1519 -273 (not in our sample) 
has an estimated core size as low as 15 to 30 //as (jMacquart et al.ll2000l ). while the most compact 
component of the extreme s cintillator J1819+3845 has been estimated to be as smal l as ~ 7uas 



(Ma cquart de Bruynll2007l ). In our present data, the rapid scintillator J1328+6221 (jKoay et al 



2011bh has the lowest upper limit of #4.9, estimated to be < 15/uas. Although Figure [T2l shows that 
the estimated lower limits of #4.9 in some sources drop well below 1/ias if their compact fractions 
are suffciently small, it is unknown if the compact fractions of these sources do indeed have values 
as low as 0.1. 



The very compact, ~ lO^uas sources give 9 sca t < 8[i&s at 4.9 GHz. Again, this can be extrapo- 
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scat ^ 3 mas at 0.33 GHz. Also, 6 scat < 264//as at 1. GHz , 



lated to 6 scat < 126^as at 1.4 GHz and 9, 
which is about a factor of 2 lower than the upper limit of 500/ias obtained by lLazio et al.1 (120081 ) 
at the same frequency. This limit is comparable to that for the sight-line to the Gamma-ray Burst 
GRB 970508 inferred from its angular size of < 3//as at 8.4 GHz ( < 9/ias at 4.9 GH z), determined 
from observations of diffractive scintillation in its radio afterglow (IFrail et al.lll997l ). 



A common parameter used to quantify the level of turbulence in the ISM is the spectral 
coefficient, C 2 , for a truncated power law distribution of electron density fluctuations (5n e ) in the 
ISM: 



2vr 



2?r 



(13) 



q is the wavenumber, Iq and l\ are the outer and inner scales of 5n e , and (3 is usually assumed 
to be 11/3 for a Kolmogorov spectrum. The scattering measure (SM), which can be derived from 
observables, is then the line-of-sight path integral of C 2 to the source at distance Dg- 



SM 



dsCl. 



(14) 



Maximum yalues of the SM f o r the IGM can be computed from the upper limits of 9 sca t, using 
the equation in iTaylor k, Corded (119931 ) . extended to cosmological scales (Koay k, Macquart, in 
prep): 



SM< 



' scatlmax) 

( Ds \ ( v \ 22 / 1 
128 mas I J VI GHz/ { ZL 



>1.2 



kpcm" 
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3 



(15) 



where in the cosmological context, Ds is the angular diameter distance to the source and Dls is 
the angular diameter distance from the source to the scattering screen in the IGM. A Kolmogorov 
spectrum is assumed for the electron density fluctuations in the IGM, modeled as a thin screen 
located at an effective distance equivalent to a redshift of zl- It is likely that scatt er broadening 



will b e dominated by nearby screens, due to the geometrical 'lever arm' effect (see iRickett et al. 



(|2007l ) for more details), as well as the redshift dependence of the rest-frame frequency at the 
scattering screen for a fixed observing frequency. This may explain why we did not observe a 
significant redshift dependence of Rd, even though the mean baryonic density of the Universe 
scales with (1 + z) 3 , and the probability of sight-lines intersecting potential scattering regions (i.e. 
galaxies and their progenitors, the intracluster medium, void walls) increases with source redshift. 
At z L ~ 0, we obtain SM < 3.3 x 10~ 5 kpcm -20 / 3 for 9 scat ( max ) ~ 8//as. 

We can also place quantit ative limits on the level of turbulence in the IGM. The SM can be 
expressed as (|Lazio et al.ll2008l ): 

SM = CsmF^Ds, (16) 
where the cons tant Csm = 1.8m~ 20 / 3 cm 6 , n e is the electron density and F is a fluctuation param- 



eter, given by (jTavlor h Corded 1 19931 ): 



rj \ lpc 



(17) 



£ is the normalized intercloud variance of the mean electron densities of each cloud, e is the nor- 
malized variance of the electron densities within a single scattering cloud, t] is the filling factor for 
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ionized clouds in the path, and In is the outer scale of the density fluctuations with Kolmogorov 
turbulence. Following Lazio et al. (j2008l ). we adopt < 2.2 x 10 _7 cm~ 3 at z ~ based on the 
assumption that helium is fully ionized, and £ ~ e ~ rj ~ 1. We therefore obtain F < 230 and 
lo > 3 x 10~ 4 pc for a source at z ~ 1 (Dg ~ 1.7 Gpc). This is consistent with what we know, 
where Iq can range from ~ 1 pc if the scattering occurs at an intervening spiral galaxy similar to 
our own, up to the ~ 0. 1 Mpc scales of t he lar gest jet sources that can inject turbulence into the 
IGM. We also note that ICordes &: Laziol (|2003l ) give F ~ 0.2 and F ~ 10 respectively in the thick 
disk and spiral arm components of our Galaxy for comparisons. 



4. Conclusions 



We analyzed data from a VLA survey of ISS in 128 sources at 4.9 GHz and 8.4 GHz to determine 
the origin of the redshift dependence of AGN ISS. We made use of tw o ISS models to interpret the 
data, one an analytical approximation ( Narayan 19921: Walkerl 19981) for the as ymptotically weak 
and strong ISS regimes, the other a fitting function ([Goodman &: Naravanll2006l ) that is applicable 
at the transition between the weak and strong ISS regimes. We also took into consideration the 
various selection effects in the source sample. We can summarize our findings as follows: 



1. The examination of the ratio of the SFs for each individual source is a good strategy for 
mitigating source selection effects in the sample, negating the redshift dependence of source 
luminosities, compact fractions, Doppler boosting factors and rest frame emission frequencies. 
Three effects can increase the ratio of -Ds.4(4d) to -D49 (4d) from Rd ~ 0.4 in the weak ISS 
regime to Rd > 1: (1) scatter broadening, (2) scintillation at the strong ISS regime, or at 
the transition between weak and strong ISS, (3) and sufficiently large scintillation timescales 
so that the SFs do not saturate at one or both frequencies so that Z?s.4(4d) rises faster than 
D 4 . 9 (4d). 

2. The examination of the correlation of the SF ratios, Rd, with source mean flux densities, 
line-of-sight Ha intensities and source redshifts allow these three competing causes of large 
Rd to be discriminated. 

3. We observed no significant scatter broadening in our sources at the scales of tens and hundreds 
of /ias probed by our survey, due either to the ISM or the IGM. We found no significant increase 
of IGM scatter broadening in the z > 2 sources compared to that of the z < 2 sources, ruling 
it out as the cause of the redshift dependence of ISS. In performing the analysis of Rd for 
the > 3cr variable sources, we note that we are including only the most variable sources at 
any redshift, which could mean that they are the least scatter broadened sources. Similar 
observations with higher sensitivity instruments such as the planned Square Kilometre Array 
(SKA), will enable Rd to be accurately estimated for sources with even lower D(t) at both 
frequencies, to determine if the sources with lower variability amplitudes are scatter broadened. 
Another weakness of the present analyses is that ~ 85% of the —0.4 < af'g < 0.4, high-redshift 
sources lie between 2 < z < 3. There is a dearth of sources at z > 3. The inclusion of more 
z > 3 sources in similar future surveys will more robustly determine if there is significant 
scatter broadening of sources beyond z ~ 3. 
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4. We infer that angular broadening in the IGM at 4.9 GHz is < HOfias for all lines of sight 
to our sources, and down to < 8/t/as for sight-lines to the ~ Wfias sources. We also obtain 
an upper limit to the scattering measure (SM) of the IGM at 3.3 x 10 -5 kpcm~ 20 / 3 for these 
latter lines of sight. 

5. We found a statistically significant steepening of source spectral indices (af'g) with source red- 
shift, which partially accounts for the redshift dependence of AGN ISS. This z-afg correlation 
can be attributed to selection effects or frustrated AGN jets in high-redshift environments. 
Follow-up high-resolution imaging of these sources using VLBI or space VLBI may help to 
discriminate between these two effects. 

6. Selecting sources in the spectral index range of —0.4 < af g < 0.4, the redshift dependence of 
AGN ISS is found to be still significant, and can be successfully modeled by a (1 + z) - 5 scaling 
of intrinsic angular sizes of a flux and brightness temperature-limited sample of sources due 
to the space-time metric of an expanding Universe. 
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A. Estimation of Transition Frequencies and Scattering Screen Distances 



This section presents the relations used to estimate the transition frequency (vt) between weak 
and strong ISS, as well as the scattering screen distance for each source. These values were used for 
the Monte Carlo simulations in Section 13.31 and to obtain estimates of the apparent angular sizes of 
the sources in Section I 



The emission measure (EM) is the integral of the square of the electron density along the path 
from the observer to the sour ce, and for the ISM is related to the line-of-sight Galactic Ha intensity 
(I a ) in units of Rayleighs as (H affner et al.lll998l ): 



EM 



2.75 T 4 °- 9 / Q cnT 



pc, 



(Al) 



where T4 is the t emperature of the io nized cloud in units of 10 4 K, typically ~ 8000 K for the warm 
ionized medium (Haffner et al.lll9 98). The transition frequency between weak and strong ISS, is 
then given bv lCordes Lazid (|2003l ) as: 



vt 



318 SM^(Asm)^ GHz, 



(A2) 



where -Dism is the effective distance to the ISM scattering screen in units of kpc, while SM is the 
scattering measure of the ISM, which is defined as the path integral of the strength of turbulence 
in the ISM along the li ne-of-sight to the sourc e (see Section [3.51 for more details on the SM), and 
has units of kpcm -20 / 3 . ICordes &: Lazid ( 20031 ) also give the relation between the SM and the EM, 
which for a thin screen can be estimated as: 



EM = 544.6 1 



2/3 
( 



! (1 + e 2 )SM pccnT 



(A3) 



where e is the normalized variance of the electron densities within the scattering cloud, which we 
assume to be ^ 1. Iq is the outer scale of the turbulence in units of pc, which has been estimated to 



This preprint was prepared with the AAS ]>Tj5X macros v5.2. 
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be < 100 pc (jHaverkorn et al.ll2008l ). We use 100 pc for our calculations. Combining Equations IA1I 
to IA31 we obtain: 



v t = 318 (Asm)" 



T9 



.9^2 



198 R 



Jo /3 (l + ^ 2 ), 



17 



GHz. 



(A4) 



The distance to the scattering screen, used in Equation IA3I and for the iGoodman Sz Narayan 

(j2006l ) fitting functions for ISS, is calculated for each source as -Dism = 0.35 x csc|6| kpc, where b is 
the Galactic latitude of the source. The value of 0.35 was selected as the constant of proportionality 
so that .Dism ~ 0.5 kpc at b = ±45°, increasing to ~ 2.0 kpc at b = ±10° and decreasing to ~ 0.35 
kpc at b = ±90°. 



It has to be noted that Equation I A3I assu mes a particular outer sc ale of Kolmogorov turbulence, 
and gives only the upper bounds of EM (see ICordes fc Lazid (|2003l )). Therefore, Equation IA4I in 
fact gives only the upper limit to i>t, and initial calculations give 4 < f£ < 80 GHz with a median 
of ~ 12 GHz. We know from Figure [5] that this is not the case, so a factor of 0.25 is multiplied to 
Uf to reduce its values to a range of 1 < v% < 20 GHz with a median of ~ 3 GHz. 
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Fig. 1. — The top panel shows a scatter plot of source spectral indices against source redshift for 
all 128 sources. The horizontal dashed lines indicate a\\ = —0.4 and a\\ = 0.4. The bottom 
panel shows the mean spectral indices in four redshift bins for all 128 sources as well as for the 
-0.4 < a 8 4 i < 0.4 sources. The error bars in the binned plots indicate la errors in the mean. 
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Table 1. List of sources and their observed properties. 



Source 


Sa q (Jy) 


Sr a (Jy) 


4.9 


Da q(4d) 


Do 4 (4d) 


R D 


ID 


z 


la (R) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


J0009+1513 


0.15 


0.12 


-0.44 


1.89e-4 ± 1.55e-4 


1.44e-4 ± 1.09e-4 


- 


fsrq 


2.2 


0.7 


J0017+5312 


0.59 


0.64 


0.14 


4.68e-4 ± 1.57e-4 


1.53e-4 ± 4.86e-5 


0.33 ± 0.15 


fsrq 


2.6 


13.3 


J0017+8135 


1.36 


1.26 


-0.13 


1.89e-5 ± 1.52e-5 


4.89e-6 ± 5.48e-9 




fsrq 


3.4 


2.2 


J0056+1625 


0.19 


0.23 


0.33 


2.40e-3 ± 8.13e-4 


1.91e-3 ± 7.00e-4 


0.80 ± 0.40 


bile 


0.2 


0.8 


J0108+0135 


1.53 


2.06 


0.56 


6.49e-4 ± 1.87e-4 


7.94e-5 ± 7.61e-5 




fsrq 


2.1 


0.7 


J0122+0310 


0.11 


0.11 


-0.04 


3.19e-4 ± 1.47e-4 


3.70e-4 ± 1.47e-4 


1.16 ± 0.70 


fsrq 


4.0 


0.5 


J0122+2502 


0.75 


0.66 


-0.21 


6.65e-5 ± 4.59e-4 


8.03e-7 ± 2.57e-4 




fsrq 


2.0 


0.9 


J0126+2559 


0.81 


0.66 


-0.39 


6.85e-5 ± 4.63e-5 


3.25e-6 ± 1.02e-4 




fsrq 


2.1 


1.0 


J0135+2158 


0.18 


0.14 


-0.37 


1.02e-3 ± 3.74e-4 


2.82c-4 ± 1.78c-4 


0.28 ± 0.20 


fsrq 


3.4 


0.9 


J0154+4743 


0.50 


0.60 


0.35 


7.53e-4 ± 1.58e-4 


9.47c-4 ± 2.09c-4 


1.26 ± 0.38 


fsrq 


1.0 


8.6 



Note. — This table is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown 
here for guidance regarding its form and content. (1) J2000.0 IAU name; (2) 4.9 GHz mean flux density; (3) 8.4 GHz mean 
flux density, may slightly differ from values published in Paper I following the removal of data points on days in which the 
4.9 GHz subarray encountered data losses, to avoid biases; (4) Spectral index, may differ from values published in Paper I; (5) 
4.9 GHz SF at time-lag of 4 days, the errors are 95% confidence bounds in fitting the model; (6) 8.4 GHz SF at time-lag of 
4 days, may differ from values published in Paper I, the errors are 95% confidence bounds in fitting the model; (7) Ratio of 
Dg.4(4d) to £>4.g(4d), given only for > 3cr variable sources; (8) Optical identification, flat spectrum radio-loud quasar (fsrq), 
BL Lac object (bile), Seyfert 1 galaxy (syfl), narrow-line radio galaxy (nlrg) or no ID available (null), obtained from the NASA 
Extragalactic Database (NED), SIMBAD database and Pursimo et al. (submitted); (9) redshift, obtained from the NASA 
Extragalact ic Database ( NED), SIMBAD database and Pursimo et al. (submitted); (10) line-of-sight Ha intensity, obtained 
from lHaffner et alj J2003h 



Table 2. Comparison of fitted Rd values at low and high redshift. 



AGN Class 


Mean Flux Density 


2 < 2 


z > 2 


Type 1 


S4.9 < 0.3 Jy 


0.39 ± 0.20 (27) 


0.70 ± 0.46 (13) 


Type 1 


S4.9 > 0.3 Jy 


0.64 ± 0.51 (10) 


1.25 ± 0.62 (11) 


Type & Type 1 


54.9 < 0.3 Jy 


0.38 ± 0.19 (29) 


0.70 ± 0.46 (13) 


Type & Type 1 


S'4.9 > 0.3 Jy 


0.83 ± 0.45 (19) 


1.25 ± 0.62 (11) 



Note. — The numbers in brackets in the third and fourth columns indicate the 
number of sources in each category. 



Table 3. Input parameters for ISS model in Section I3.3L 



Parameter 




Symbol 


Value 


Scattering screen distance from Earth 




-DlSM 


500 pc 


Scattering screen velocity 




v s 


50 kins" 1 


Transition frequency between weak and strong 


ISS 


Vt 


4.0 GHz 


Source intrinsic brightness temperature 




Tb,int 


10 11 K 


Source compact fraction 




fc 


0.5 


4.9 GHz and 8.4 GHz mean flux density (stroni 


5 sources) 


s v 


1.00 Jy 


4.9 GHz and 8.4 GHz mean flux density (weak 


sources) 


S v 


0.15 Jy 
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iog 10 [D^(4d)] '°g 10 [ D 8.4( 4d >] 



Fig. 2. — Distribution of Z?4.g(4d) (left) and Ds.4(4d) (right) in the low (top) and high (bottom) 
redshift sample of sources. Only the 102 -0.4 < af; 4 . < 0.4 sources are shown, classified into Type 
or Type 1 AGNs based on their optical IDs. Since D no i se « 1 x 10~ 4 on average at both frequencies, 
we include all sources with D(4d) < 1 x 10" 4 in the -4 < logi [-D(4d)] < -3.5 bin. The high 
redshift sources are significantly less variable than their low redshift counterparts. 
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Fig. 3. — 4.9 GHz spectral luminosities plotted against source redshifts, with sources separated into 
four Z?4.g(4d) bins. Only the 102 —0.4 < af g < 0.4 sources are included. The solid curve gives 
the L4.9-Z relation for a 0.3 Jy source, separating the strong (S4.9 > 0.3 Jy) sample of sources from 
the weak (54.9 < 0.3 Jy) sample of sources. The vertical dashed line indicates the redshift cutoff of 
Z?4.9(4d) at z ~ 2, while the dashed horizontal lines indicate possible luminosity cutoffs of Z?4.g(4d) 
at L4.9 ~ 10 28 WHz" 1 and L4.9 ~ 10 27 WHz~ 1 for the strong and weak sources respectively. 
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Fi g. 4. — The structure f u nction ratio, Rd, calculated for j3 = 0, 1 and 2.2 using the fitting function 



in 



Goodman fc Narayanl (|2006l ) and plotted against the transition frequency between weak and 
strong scattering. We have adopted values of v s = 50kms , 6 = 100 ^as at 4.9 GHz, and a 
scattering screen distance of -Dism = 500pc. The vertical dashed lines represent v = 4.9 GHz and 
v = 8.4 GHz respectively. 




Fig. 5. — Histogram of Rd for 72 sources with D(4d) > 3c at both frequencies, classified as Type 
or Type 1 AGNs. The dashed vertical lines denote Rd = 0.4 and Rd = 1.8. 
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Fig. 6. — Plot of Z?8,4(4d) against Z?4.g(4d) in logarithmic (left) and linear (right) scales with sources 
classified as 'weak' or 'strong' based on their observed mean flux densities at 4.9 GHz, S4.9. Only 
sources with D(4d) > 3a at both frequencies are included. The dotted lines represent Rd = 0.4 
and Rd = 1-8. The solid lines represent linear fits to the #4.9 < 0.3 Jy and 54.9 > 0.3 samples, 
while the dashed lines represent 99% confidence bounds for those fits. 
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Fig. 7. — Plot of Z?8.4(4d) against D4.g(4d) in logarithmic (left) and linear (right) scales, where 
the sources are classified based on their line-of-sight Ha intensities, I a , in units of Rayleighs. Only 
sources with D(4d) > 3cr at both frequencies are included. The dotted lines represent Rd = 0.4 
and Rd = 1.8. The solid lines represent linear fits to the I a < 0.5 R and I a > 0.5 R samples, while 
the dashed lines represent 99% confidence bounds for those fits. 
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Fig. 8. — Plot of 1^8. 4 (4d) against L>4.g(4d) in logarithmic (left) and linear (right) scales, where the 
sources are classified based on their redshifts. Only sources with D(4d) > 3a at both frequencies 
are included. The dotted lines represent Rd = 0.4 and Ro = 1-8. The solid lines represent linear 
fits to the z < 2 and z > 2 samples, while the dashed lines represent 99% confidence bounds for 
those fits. 
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Fig. 9. — Observed D(4d) at 8.4 GHz (left) and 4.9 GHz (right) plotted against redshift, shown 
as scatter plots (top) and in bin averages (bottom), for both the weak (S4.9 < 0.3 Jy) and strong 
(1S4.9 > 0.3 Jy) sources. The vertical error bars in the binned plots represent one standard error in 
the mean. The lines in all panels show computed model values of D(4d) for various values of the 
source Doppler boosting factor, assuming that cosmological expansion leads to a (1 + z) 5 scaling 
of the intrinsic angular diameter in a flux and brightness temperature-limited sample of sources. 
The model parameters used are listed in Table El 
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Fig. 10. — Observational values of the fitted Rd in two redshift bins with their corresponding o»7o 
confidence bounds la errors), separated into weak (54.9 < 0.3 Jy) and strong (S4.9 > 0.3 Jy) 
sources. These are shown together with model values of Rd for various source Doppler boosting 
factors. For the blue and red curves, the model parameters are given in Table [3l The black curves 
show the corresponding model values for the strong sources at a screen velocity of 20 kms -1 with 
the other parameters unchanged. 
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Fig. 11. — Observed values of -D4.g(4d) (left panel), _D4.g(4d) (middle panel) and Rd (right panel) 
plotted against their respect ive model values obtained by applying the Monte Carlo method to the 



Goodman Naravanl (|2006l ) fitting functions (described in Section [373]) . For D^_g(4d) and Dg.4(4d), 



all 102 —0.4 < CK49 < 0.4 sources are shown. For Rd, only the 72 sources with > 3a variability are 
shown. The solid diagonal lines show where the observed values of -D(4d) and Rd are equal to their 
model values. 
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Fig. 12. — Top panel: Constraints on the 4.9 GHz apparent angular sizes, #4. 9, for all 128 
sou rces (including those with a fj > 0.4 and afg < —0.4), calculated based on Z?4.g(4d) using 
the lGoodman &; Naravanl (|2006l ) fitting function (see Section [3.51 for more details). Bottom panel: 
Upper limits of #4.9 for all sources in which they can be obtained, shown proportional to the sizes 
of the circles and plotted in Galactic coordinates. They are also effectively upper limits of 9 sca t for 
all lines of sight to our sources. The circles are color-coded based on the redshifts of the sources. 
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Fig. 13. — Ru at increasing values of SC at for sources with intrinsic angular sizes of 10 fias and 100 
fias. The plots are based either on model calcu lations in the weak ISS regim e assuming that the 
SFs have all saturated, or the fitting formula of iGoodman fc Narayanl ( 200a ) with v s = 50kms _1 , 
Dism = 500pc and v t = 4.0 GHz. The dotted horizontal line represents Rd = 0.57, which is the 
upper limit to Rd fitted to all sources in the low redshift sample, at a confidence level of 99%. 
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